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Abstract 

At large N c the nucleon can be viewed as a soliton of the effective chiral lagrangian. 
This picture of nucleons allows a consistent nonperturbative calculation of the leading- 
twist parton distributions at a low normalization point. We derive general formulae 
for the transverse spin quark distribution hi (x) in the chiral quark-soliton model. We 
present numerical estimates and compare them to the results obtained in other models. 
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1. The transverse polarization distribution of quarks introduced by Ralston and Soper 
[[J and usually called h\{x) has gained an increasing amount of attention of theorists [0. [|. 
H> i; H 0- The perspective of the experimental measurement of h\(x) || |j has stimulated 
attempts to calculate hi(x) in various models of strong interactions 0, |5], 10, [XT . 



Recently a new approach to the calculation of quark distribution functions has been 
developed [O within the context of the effective chiral quark-soliton model of nucleon |T3 



In this paper we apply this approach to the calculation of the transverse spin distribution 
function hi(x). Although in reality the number of colours N c = 3, the academic limit of large 
N c is known to be a useful guideline. At large N c the nucleon is heavy and can be viewed as 
a classical soliton of the pion field [14, 15fl. In this paper we work with the effective chiral 



action given by the functional integral over quarks in the background pion field [16|, [17], [18 

exp (iS eS [ir(x)\) = f DipD^ exp (i f d 4 x^{i$ - MIP 6 )^ 



U = exp {m a (x)T a } , W 5 = exp [m a {x)T a lh \ = ^~^U + — ^£7+. (1) 

Here ip is the quark field, M is the effective quark mass which is due to the spontaneous 
breakdown of chiral symmetry and U is the SU(2) chiral pion field. The effective chiral 
action given by eq. ([!]) is known to contain automatically the Wess-Zumino term and the 
four-derivative Gasser-Leutwyler terms, with correct coefficients. Eq. ([!]) has been derived 



from the instanton model of the QCD vacuum [Tj| [19), which provides a natural mechanism 
of chiral symmetry breaking and enables one to express the dynamical mass M and the 
ultraviolet cutoff intrinsic in eq. (H) through the Aqcd parameter. It should be mentioned 
that eq. (|1|) is of a general nature: one need not believe in instantons and still can use eq. (|1|). 

An immediate application of the effective chiral theory ([[]) is the quark-soliton model 
of baryons of ref. |Tj|, which is in the spirit of the earlier works |2^, ^l[| . According to this 
model nucleons can be viewed as iV c (=3) "valence" quarks bound by a self-consistent pion 
field (the "soliton" ) whose energy coincides in fact with the aggregate energy of the quarks of 
the negative-energy Dirac continuum. Similarly to the Skyrme model large N c is needed as a 
parameter to justify the use of the mean- field approximation, however, the 1/N C corrections 
can be and in some cases have been computed [f22], |24], |26|, |27J . The quark-soliton model 



of nucleons developed in ref. jn| includes a collective-quantization procedure to deal with 
the rotational excitations of the quark-pion soliton. 

Turning to the calculation of the transverse spin distribution function hi(x) we note that 
the model possesses all features needed for a successful description of the nucleon parton 
structure: it is an essentially quantum field-theoretical relativistic model with explicit quark 
degrees of freedom, which allows an unambiguous identification of quark as well as antiquark 
distributions in the nucleon. This should be contrasted to the Skyrme model where it is not 
too clear how to define quark and antiquark distributions. 

2. The transverse spin quark distribution function for flavour / is defined as follows [5| 

h f M = J —e iX *(P,S T $ f m l5 $ T M*n)\P,S T ) . (2) 
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Below we work in the nucleon rest frame where the transverse nucleon spin St and the 
light-cone vector n can be chosen as follows: 



rv 



1 



M, 



% -1,0,0), ^ = (0,0, 0,1). (3) 



N 



The operators in eq. (|2|) depend on the QCD normalization point /i. In contrast to QCD, 
the effective chiral field theory is nonrenormalizable and contains an explicit ultraviolet cut- 
off. In the instanton vacuum model this cutoff appears as the inverse average instanton size: 



p 600 MeV [0, |25[ . The results obtained below thus refer to a low point normalization 



point p of the order of 600 MeV. 

Now we can rewrite (0) in the form 



h[( x ) = — [ dz°e ixMNZ ° 
An J 

— oo 

x (P, Sr|l#(0)(l + TVhsT^/WlP, St>[ 1= _^ 2 _3 =0 • (4) 

Let us remind the reader how the nucleon is described in the effective low-energy theory 
(|Tj). Integrating out the quarks in (|I|) one finds the effective chiral action, 

S eS fr a (x)} = -N c Sp\ogD(U), D{U) = id -H(U), (5) 

where H(U) is the one-particle Dirac hamiltonian, 

H(U) = -«7 V<9 fc + M 7 °[/ 75 , (6) 

and Sp . . . denotes a functional trace. 

For a given time- independent pion field U = exp(z7r a (x)r a ) one can determine the spec- 
trum of the Dirac hamiltonian, 

H$ n = EA. (7) 

It contains the upper and lower Dirac continua (distorted by the presence of the external 
pion field), and, in principle, also discrete bound-state level(s), if the pion field is strong 
enough. If the pion field has unity winding number, there is exactly one bound-state level 
which travels all the way from the upper to the lower Dirac continuum as one increases the 
spatial size of the pion field from zero to infinity [I3|. We denote the energy of the discrete 



level as Ey ev , —M < Ey ev < M. One has to occupy this level to get a non-zero baryon 
number state. Since the pion field is colour blind, one can put iV c quarks on that level in 
the antisymmetric state in colour. 

The limit of large N c allows us to use the mean-field approximation to find the nucleon 
mass. To get the nucleon mass one has to add N c E\ ev and the energy of the pion field. 
Since the effective chiral lagrangian is given by the determinant (|5]) the energy of the pion 
field coincides exactly with the aggregate energy of the lower Dirac continuum, the free 
continuum subtracted. The self-consistent pion field is thus found from the minimization of 



the functional 
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M N = min N c E leY [U] + £ (E n [U\ - E<®) . (8) 

[ E n <0 J 

From symmetry considerations one looks for the minimum in a hedgehog ansatz: 

U c (x) = exp [z7r a (x)r a ] = exp [in a r a P{r)} , r = |x|, n = -, (9) 

where P(r) is called the profile of the soliton. 

The minimum of the energy is degenerate with respect to translations of the soliton 
in space and to rotations of the soliton field in ordinary and isospin space. For the hedgehog 
field @ the two rotations are equivalent. Quantizing slow rotations of the saddle-point 
pion field [15, 13] leads to the projection on a nucleon with given spin (£3) and isospin (T 3 ) 



components which includes the integration over all spin-isospin SU(2) orientation matrices 
R of the soliton, 

(S = T, S 3 , T 3 \ . . . \S = T, S 3 , T 3 > = JdR ^ T (R) . . . </> s s ^(R) • (10) 

Here 0f^(i?) is the rotational wave function of the nucleon given by the Wigner finite- 
rotation matrix [|H| 13|[ : 



s t(R) = V2STT(-lf+ T W s _^ S3 (R). (11) 



J S 3 T 3 \ 



3. Eq. (0) contains a matrix element of a non-local quark bilinear operator in the nucleon 
state with definite 4-momentum P and spin and isospin components. According to [Q one 
can write a general expression for such matrix elements; the time dependence of the quark 
operators is accounted for by the energy exponents. Taking nucleon at rest, we can write 
this matrix element as a sum over all occupied quark states 

(S 3 ,T 3 \i]} t (x°,x)r j (y°,y) \S 3 ,T 3 ) = 2M N N C J d 3 X J dR<P* SsT3 (R) 

x J2 exp[zE n (x° - y )] ^(x - X) (R^' R 9 g> $^'(y - X) fa^R) . (12) 

n 

occup. 

Here we have written explicitly all the flavour (/, g = 1,2) and the Dirac = 1,...,4) 
indices for clarity. The functions $ n are eigenstates of energy E n of the Dirac hamiltonian 
(^]) in the external (self-consistent) pion field U c . Summation over colour indices is implied 
in the quark bilinears, hence the factor iV c . 

Applying the general formula of the effective chiral theory ([T2] ) to the calculation of the 
transverse spin distribution function (^) we obtain in the leading order of the 1/N C expansion 

00 

N r M, 



2,71 J J n 

~°° occup. 

x(n\R^T f R(l + 7 V)757 3 exp(-^V)|n). (13) 
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Here T* is the matrix projecting onto the flavour /. We use notation p k = —idk for the 
operator of the 3-momentum. One can easily check that in the leading order of the 1/N C 
expansion only the isovector polarized distribution survives: 



2n J J „ 



■iE n z° 



occup. 



x (n\RW 3 R(l + 7V)757 3 exp(-i*V) \n) . (14) 
The integral over the orientation matrix R can be easily computed. As a result we obtain 

oo 

N cMN_ f j^ JxM N z° ST e -iE n z° 



hUx) - hf(x) = / dz°e ixMNZ ° Yl 



n 

oo 



occup. 

x(n|r 3 (l + 7°7 1 )757 3 exp(-^V 1 )l^) • (15) 

Taking into account the hedgehog symmetry (|]) we can replace here 1 <->• 3. Passing to the 
momentum representation and integrating over z°, we obtain the expression convenient for 
practical calculations: 

N c Mn f d 2 p ± 



occup. 

x(n|p)(l + 7 7 3 )757V(p|n) . (16) 

p i =xM N —E n 

It is well known that the transverse spin distribution function hi satisfies the sum rule 

dx h{{x) = (17) 

-l 

where g T is the tensor charge of the nucleon defined as follows 

(P,S$ f [>y fl ,>y»]il>f\P,S) =g f T u(P,S)[ lfinv ]u(P,S). (18) 

Here u(P, S) is the standard Dirac spinor describing the nucleon with momentum P and 
spin S. 

Substituting expression for h\{x) (0) into the lhs of the sum rule (fPTD and extending the 
integral over x to the whole real axis (which is justified by the large N c limit [ |T2[] ) we obtain 

r N 

dxh{{x) = ^ («l757 1 T 1 |n). (19) 

J 6 n 

— -L occup. 



The result in the rhs coincides with the expression obtained in paper |26| for the tensor 
charge of nucleon so that the sum rule fll7D automatically holds in the effective chiral theory. 



4. In eq. ( |16D one has to sum over all occupied levels of the Dirac Hamiltonian, including 
both negative Dirac continuum and discrete level. 
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In contrast to the longitudinal spin quark distribution function gi(x) whose Dirac con- 
tinuum contribution is ultraviolet divergent ||12|| , in the transverse case the continuum con- 
tribution is ultraviolet finite. This can be seen from the gradient expansion for hi(x) where 
the first nonvanishing term appears in higher orders compared to the gradient expansion of 
gi(x). This also hints that the continuum contribution to hi(x) is relatively small, whereas 
the continuum contribution to the longitudinal spin distribution function g\{x) is rather 



sizable [fL2 |. An additional argument in favour of the suppression of the continuum contribu- 
tion to h\ comes from the numerical calculation of the tensor charge performed in [T2^, p7| , 
which shows that the contribution of the Dirac continuum is essentially smaller than the 
contribution of the discrete level. Since the sum rule ( |l7|) holds in our model, it is natural to 
expect that the discrete level gives the dominant contribution also to the quark distribution 
function h\{x). 

The bound-state level occurs [[L3|] in the grand spin K = and parity II = + sector of 
the Dirac hamiltonian (||). In that sector the eigenvalue equation takes the form 



M cos P(r) 



d_ 

V dr 



- M sin P(r) 



d 2 

-- M sin P{r) 

or r 

-M cos P(r) 



h {r) 



J 



lev 



ho{r) 



We assume that the radial wave functions are normalized by the condition 



(20) 



drr 2 [h 2 (r) + jf(r)\ 



1. 



(21) 



Introducing the radial wave functions in the momentum representation 



h(k) = I dr r 2 ho(r)]l — Ji(kr) , j(k) - ' 

o o 

we can write the expression for the level contribution to hi in the form 



drr^jiMil-Jafkr), (22) 



h%(x) - hf(x) 



= -N C M N 

lev 3 

\xM N -E lcv \ 



dk ^ 2{k) + (xM N - E lev ) 2 m 



k 2 



(23) 



5. Recently the following inequality has been suggested by Soffer [28 

q f (x) + g((x) > 2\h{(x)\. 



(24) 



Its status in QCD has been studied in Let us check this inequality in the chiral quark 
soliton model. A delicate point is that in the leading order of the 1/N C expansion the 
unpolarized quark distribution qf is saturated by the singlet part whereas g{ and h[ get the 
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main contribution from their nonsinglet parts. Therefore in the large N c approach the above 
inequality can be rewritten as follows: 



[q u {x) + q d (x)] - \g»(x) - g£(x)\ > 2\K}(x) - hf(x)\ . (25) 

Note that the structure of the expressions for all quark distribution function in the effective 
quark soliton model entering inequality (|25|) is essentially the same. They contain a sum of 
diagonal matrix elements over occupied states of the Dirac hamiltonian, differing only by 
the spin and isospin matrices T appearing in these matrix elements: 

[q u (x) + q d (x)], r g = (l + 7 V), 

\g?(x) - gf(x)} , r L = + 7 °7 3 )7 5 r 3 , 

[h^x)-h({x)}, r T = l(l + 7°7 3 )757V 1 (26) 



(the explicit expressions for q* (x) and g\{x) can be found in ref. |12||). Taking into account 
the following inequalities (understood in the matrix sense) 



|(l + 7 V)7 5 r1 <1 + 7V, (27) 
|(1 + 7 ° 7 3 )757V| < l+ 7 °7 3 (28) 



one can easily see that the Soffer inequality ( Ell ) holds for the contribution of each separate 
occupied level. 

One should keep in mind that in our model the quark distribution functions q(x) and 
gi(x) are ultraviolet divergent and in principle the regularization can violate the above 
argument. However, at least in the limit of large cutoff we can prove the validity of the 
Soffer inequality. Indeed, function h\(x) is ultraviolet finite whereas q(x) and gi(x) are 
ultraviolet divergent. The analysis of these divergences has been performed in ref. |12[ and 
it follows from this analysis that the ultraviolet divergence of the combination [q u (x) + 
q d (x)} — \gi(x) — gf{x)\ is given by an explicitly positive expression so that in the limit of 
large cutoff this expression definitely becomes larger than the ultraviolet finite hi(x), which 
proves that the Soffer inequality holds in this limit. 

6. We have calculated numerically the level contribution to the transverse isovector 
quark distribution. In this calculation we use the parameters obtained by using a variational 
estimate of the profile of the pion- field soliton (see eq. (§)) performed in ref. [D| yielding 
for M = 350 MeV 



r 



P(r) = -2 arctan -£ , r w 1.0/M, M N « 1170 MeV. (29) 



This profile function has a correct behaviour at small and large distances and is stable with 
respect to small perturbations. In our numerics we use the analytic profile (]29f) . 

The results of our calculations for g\ and hi are presented at Fig.l. The antiquark dis- 
tributions are shown at Fig.2. From the latter plot we can see that the antiquark transverse 
spin distribution is essentially smaller than the corresponding longitudinal distribution which 
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is related to the fact that the longitudinal distribution gets a large contribution from the 



quarks of the Dirac continuum [12| in contrast to the transverse distribution. 

Comparing our results to the calculations in other models - bag model [|, [|, QCD sum 
rules chiral chromodielectric model |TT], one can see that the main source of difference 
is the 1/N C approximation used in our calculations. As it was mentioned above, in the 
leading order of the 1/N C expansion we have h\(x) = —hf(x). This relation is much more 
general than the model we used and also holds in the large N c QCD. It is well known 
that typically the large N c counting is a good guide line for understanding nonperturbative 
physics. In order to understand the reliability of the 1/N C expansion in the case of hi it 
were instructive to compute the next-order 1/N C corrections. Unfortunately, the calculation 
of these corrections in our model is rather complicated in the case of the quark distribution 
functions. However, these 1/N C corrections have been computed for the tensor charge ]26], [27 



and this calculation shows that these corrections really enhance g% and suppress g^. This 
means the 1/N C corrections of our model can reconcile the large N c approach with the model 
calculations of § §, 0, [TTJ]. 
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2.0 




Figure 1: The transverse spin quark distribution function h\ — hf (solid line) plotted 
against the longitudinal quark distribution — gf (dashed) computed in the same model. 
The contribution of the discrete level to — gf is depicted by the dot-dashed line. 
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